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Abstract -We find a nontrivial regime of spinfoam quantum gravity that reproduces classical 
Einstein equations. This is the double scaling limit of small Immirzi parameter (gamma), large 
spins (j) with physical area (gamma times j) constant. In addition to quantum corrections in the 
Planck constant, we find new corrections in the Immirzi parameter due to the quantum discreteness 
of spacetime. The result is a strong evidence that the spinfoam covariant quantization of general 
relativity possesses the correct classical limit. 



Introduction. — Spinfoams P~H3] are a tentative co- 



^Jyariant quantization of general relativity. They pro- 
vide transition amplitudes between quantum states of 3- 
geometry, in the form of a Misner-Hawking sum over vir- 
^_tiual geometries jH|5j. In the so-called 'new models' [MS], 
qq intermediate quantum states are the ones of canonical loop 
quantum gravity, the SU(2) spin-network states, a rc- 
O^l markable feature that promotes the spinfoam framework 
^Jtjo a tentative path integral representation of loop quan- 
OOtrum gravity. 

, The physical picture emerging from the spinfoam grav- 
| ity is the following: spacetime is a quantum foam of virtual 
. . geometries with a discrete and purely combinatorial, rela- 
£^fcional structure, where the Planck scale plays the role of 
k> 4 natural minimal length. 

■ General relativity is described by the Hoist action 

a ■ 



s 



i 



(e A e) A F(w) + - / (e A e) A F{w) (1) 



with gauge group 50(1,3) (or 50(4) in the Euclidean 
signature), where the second term vanishes by using the 
equations of motion. The first term is the Einstein- 
Hilbcrt-Palatini action for general relativity in terms of the 
cotetrad e and the connection uj, regarded as independent 
variables (F is the curvature). The real number 7 7^ is 
called Barbero-Immirzi, or Immirzi parameter, and con- 
trols the spacetime discreteness as it enters the discrete 
spectra of area and volume operators [9, 10 . Spinfoam 
models provide a Feynman path integral, or state sum, 
based on a discretization of ([T]) over a 2-complex (gener- 
alized triangulation) and the full degrees of freedom are 
recovered in the infinite refinement limit or equivalently 



summing over 2-complexes [11] , The spinfoam theory is 
sufficiently simple [TJ] and possesses the correct symme- 
tries [T5] , Recently it has been successfully coupled to 
matter fields |14I15| . Furthermore, the large distance anal- 
ysis was able to extract the correct low-energy physics in 
some simple cases [TBHH]. 

A major open problem is to show that the spinfoam 
Feynman path integral is able to reproduce Einstein equa- 
tions when a semiclassical expansion is performed. What 
we would like to have is the spinfoam version of the fol- 
lowing semiclassical expansion of the gravitational path 
integral 



Dg^e 



(2) 



in the classical limit h — > 0, where Seh is the Einstein- 
Hilbcrt action for general relativity and on the right hand 
side it is evaluated on the classical solution g^ v of the 
equations of motion determined by the boundary condi- 
tions on the metric field g^ u . Here we propose a solution 
to this problem. 

In this letter we improve and generalize previous argu- 
ments |20j to a general 2-complex. In the continuous area 
spectrum limit 7 — > and in the semiclassical limit h — > 
we find the analogous of the WKB expansion @ for spin- 
foam quantum general relativity. We can state the regime 
in a more suggestive fashion by introducing two physical 
scales. One is the length scale of quantum gravity, identi- 
fied with the Planck length 



QG 



h 
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the other is the scale of loop quantum geometry, that is spin 1/2 representation and we have the following con- 
the scale where we can 'see' the discreteness of spacetime straintjj on the spin labels: 



LQG 



p- 



(4) 



Thus the regime of the spinfoam path-integral we look for 
is expressed by the following relation: 



l LQG 



(5) 



where I is the typical linear scale of each 4-simplex in 
the spacetime triangulation. The analysis is partly based 
on the path integral formulation of references [2"TII22] . 
Throughout the letter we work in natural units G = h = 
c = 1, but will restore some of those constants when 
needed. In particular, restoring only h the Planck length 
is lp — \ph. 

The spinfoam amplitude. — We consider the spin- 
foam amplitude [5][5] for a 2-complex a without matter. 
We restrict for simplicity to the Euclidean signature and 
to Barbero-Immirzi parameter < 7 < 1, where formulas 
get simpler. For each face / (the 2-cells) of the 2-complcx 
a there is an associated integer spin jf . Faces are oriented 
and bounded by a cycle of edges e (the 1-cells) . Each edge 
bounding a face has a source vertex s(e) and a target ver- 
tex t(e), where source/target is relative to the orientation 
of the face. To each edge let us associate SU(2) elements 
n e f (/ runs over the faces meeting at the edge e) and two 
source/target Spin(A) ~ SU(2) x SU{2) ~ 50(4) gauge 
group variables <? e . s ( e ), g e ,t(e)- The variables n e f can also 
be interpreted as unit vectors n e f in R 3 , up to a phase 
ambiguity, by saying that n is a rotation that brings a 
reference direction to the direction of n. 

The spinfoam amplitude, or partition function, for the 
2-complex a in the Bloch coherent state basis [33J is de- 
fined as 



Z 



E 
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dg ve 



dn ef Y[P f . 
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(6) 



The sum is over spins jf € N, and the integrals are over 
the Spin{^) gauge variables and SU{2) variables labeling 
the edges 1^1 The face amplitude Pf is given by 



P f =t r U eef P+ 



ef 



(7) 



where II denotes the ordered product (according to the 
cycle of edges) and 



P, 



fir 



s(e)l 



if (7 



e/l 
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e,t(e) 



(8) 



Here |n) is the SU(2) Bloch coherent state [24] for angular 
momentum^ along the direction of n, in the fundamental 

1 Wc dropped a face normalization factor which is usually taken 
as the dimension d(jf) = 2jf + 1 of the S'f/(2)-irreducible Hilbert 
space. Other normalizations arc possible but irrelevant in the present 
analysis. 

2 We introduced the notation [j, n) for the standard antilinear map 
t applied to \j,n). In the standard basis, it is given by the symbol 



r 



1±7 
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We have split the SO (A) variables g into selfdual and an- 
tiselfdual rotations (g + ,g~) € SU{2) x SU{2). The par- 
tition function © can be written in the form of a path 
integral for an action as 



Z =y^l dg ve / dn, 
in) 

5 = E 5 / = E lnP / 



(10) 

(11) 



The continuum 7—^0 limit. — The continuum limit 
of the theory is defined as the infinite refinement of the 2- 
complex [TT] (possibly undergoing a second order phase 
transition to a smooth spacetime). Differently, here we 
define a continuum limit that is suitable to describe gravi- 
tational physics with a truncation of the theory on a finite 
cellular structure of spacetime, namely with an approxi- 
mation of the full theory. The spinfoam amplitudes associ- 
ated to finite graphs can be viewed as effective amplitudes 
obtained after a coarse-graining procedure |25U26| applied 
to the spinfoam infinite 'lattice'. But since the truncated 
amplitudes are not fundamental, there is no reason to keep 
fixed the Immirzi parameter to its 'bare' value 70. The 
possibility of a renormalization of the Barbero-Immirzi pa- 
rameter has been recently advocated in different contexts 

EHES]. 

Here we explore the possibility of a running towards 
zero, 7 — > 0, simultaneously with the large-spin regime of 
the theory. Thus we consider 



j — > 00, 7—>0, j"f = const. 



(12) 



and 7j is the macroscopic physical area in Planck area 
units. Notice that the Immirzi parameter controls the 
spacetime discreteness. In particular it controls the area 
gap and the spacing between area eigenvalues, thus the 
limit (|12[) is the continuum limit for the area operator. 

In order to see the effect of (Tj"2"]) on the partition func- 
tion, let us restrict our attention to a 2-complex which is 
dual to a simplicial triangulation: vertices are dual to 4- 
simplices, and edges are dual to tetrahedra. The analysis 
parallels the one of [221129] at fixed, large spins. We are in- 
terested in making explicit the dependence on the Immirzi 
parameter, so let us decompose the action pT]) using © 
in the following way 

S f = a f (In Pf - In Pj) + -af (In Pf + In PJ ) (13) 

where we have defined the area of the triangle dual to / as 
a f — 13 f j an d Pf 1 is the face amplitude in the fundamental 



m. — m 



3 Because of equation (J9J 7 is quantized to be rational. This 
restriction is not present in the Lorcntzian version of the theory. 
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representation. We want to evaluate the partition function 
(110D in a region of macroscopic areas a / and in the limit 
7 — > with a/ fixed (so jf — > oo). Collecting all the face 
terms (1131). let us write the full action as 



s = s° + -S'. 

7 



(14) 



In the limit 7 — > with a / fixed the partition function can 
be approximated with an integral over continuous areas 



Z 



daf dg ve / dn e fe 1 



(15) 



For the stationary phase evaluation of (fT5"]) we have to take 
variations of the second term S" proportional to the large 
parameter - — > 00 . The action S' is complex with nonpos- 
itive real part Re S" < 0, so the main contribution to the 
integral comes from the critical points, where ReS" = 0. 



is an angle with the interpretation of a torsion degree 
of freedom. But the on-shell discrete connection g^ v is 
torsion-free (it is the discrete spin-connection) and this 
angle vanishes (see also [22] for more details on torsion). 

The critical points are in one-to-one correspondence 
with 4-dimensional spacetime triangulations (Regge man- 
ifolds) , where the area^l of the triangles are specified by 
a f , and their 3d normals by n e f . The Regge manifolds are 
endowed with a continuous, piecewise flat metric where 
curvature is distributional and concentrated on triangles. 
From a given critical point, one can reconstruct the area 
bivectors A v f of the triangles dual to /, in the frame of the 
4-simpliccs v, namely the metric of the Regge manifold. 

In the 7 — > expansion the action must be evaluated at 
a critical point, and using C|18|) we have 



This condition holds for ReP^" 



S\cnt. = iS R = i 2J a/O/ 
/ 



Re P~ 



0. One can 



easily show that the critical points are the solutions to where we have defined the deficit angle as 



9tv n ef = 
9ev™ef 



— 



(16) 
(17) 



where e, e' are adjacent edges in the face /, sharing the 
vertex v. If there are no solutions, the amplitude is expo- 
nentially suppressed. Now using (fT6"|). (fTT|) we have that 
the brackets 



( n ef\(gte) 1 9ve>l n e'f) 



if? 3 



(18) 



reduce to simple phases, on the critical points (on-shell). 

Furthermore, we must require that the critical points 
are also stationary. Varying S' with respect to 5*0(4) 
group variables, and evaluating at a critical point we get 
the condition 



,s'U. = o 



E 

fee 



ajUef = 



(19) 



which expresses the closure relation for the tetrahedron 
dual to the edge e. Variation with respect to the unit 
vectors n e f does not give further information, because it is 
automatically satisfied. Finally, we have to take variations 
of 5" with respect to the areas a/, but this docs not give 
further restrictions since 



dS' 
da f 







(20) 



is automatically satisfied on the critical points. Indeed 
using ([15]) one can show that 



dS' 
da f 



crit. v< zf 



K f ) = ©/ 



(21) 



4 We are only interested in the oscillatory behaviour of the inte- 
gral, so we drop one global factor I/7 per each face, that takes into 
account the measure of the Ricmann sum. 



0/ 



vef 



(22) 



(23) 



The deficit angle O / is the spacetime curvature concen- 
trated on the triangle dual to the face /. The equation 
(]22|) is the Regge form [3"3"ll3"l] of the action £*# for general 
relativity. The asymptotic approximation of the integral 
(115p requires to sum over all critical points. In general 
this will take the form of a continuous sum over the set 
C of critical points (critical manifold) , with some measure 
[i that can be computed by standard asymptotic analysis 
tools. The result is 



Z 



d[i(a f ,n ef ,g ev )e iSR +0(7) 



(24) 



where 0(7) denotes the 7-corrections to the partition 
function coming from the next orders of the asymptotic 
approximation. 

The semiclassical h — > limit. — More interest- 
ingly, let us parametrize the previous integration (1241) 
using length variables. Given that the critical points 
(af,n e f,g ev ) correspond to Regge triangulations, there 
exists an assignment of lengths l s to the sides s of the 
triangles such that the areas a/ coincide with the areas 
computed out of the lengths: a/ = a/(Z s ). Then we can 
parametrize the critical manifold C with the set of side 
lengths l s . Restoring the h dependence, let us rewrite 
(E§ as 



Z = J dfi(l s )ei s ^ +0( 7 ), 



(25) 



5 It is well known that a generic assignment of areas does not 
correspond to any Regge triangulation I3UH32I . In other words, in 
general there exists no assignment of lengths l B for the sides s of 
triangles such that aj = aj(l B ). The critical equations select only 
the Reggc-likc area configurations. 
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and the Regge action is now explicitly a function of the 
lengths 

S R (l s ) = J2a f (l s )@ f (l s ) (26) 
/ 

as in the original classical formulation [33] ■ The last ex- 
pressions are a good starting point for taking the semi- 
classical limit. Before discussing this, let us make a few 
comments on the effect of sending the Immirzi parameter 
to zero. The remarkable consequence of the continuum 
limit 7 — > we have performed in the last section is that 
the spinfoam amplitude reduces effectively, that is up to 
7-corrections, to a quantization of Regge gravity [35] given 
by formula (|25[) where the fundamental variables are con- 
tinuous lengths. The result resonates with the recent find- 
ings in the computation of the graviton propagator within 
loop quantum gravity |16H18j . In particular, as shown in 
[18] , the leading order (in the h expansion) graviton prop- 
agator G(x, y) presents the same kind of 7-corrections, 

G(x, y) = R+ . lX + f Y + h-corv. (27) 
\x - yY 

and only in the limit 7 — > the tensorial structure of the 
2-point function matches with the matrix of correlations 
R computed in quantum Regge gravity, and, even more 
interestingly, with the one given by standard perturbative 
gravity on flat space. In retrospect, at the light of the 
present general analysis the previous result (|27[) is much 
more clear. Moreover, similar 7-corrections were found by 
Bojowald in the cosmological context [36) . 

Suppose now we are interested in the semiclassical ex- 
pansion of the spinfoam amplitude (|2~5j) . This corresponds 
to looking at areas which are macroscopic, that is large as 
compared to the Planck area, 

Vr » 1, (28) 

Lp 

or equivalently to the standard WKB expansion h — > 0. 
This regime can be selected by appropriate semiclassical 
boundary conditions in the transition amplitudes, as ex- 
plained in the next section. 

The classical equations of motion are obtained by vary- 
ing (|21)|) with respect to the lengths. Using also the Schlafli 
identity [33] , which tells that the variation of the deficit 
angles do not contribute to the total variation of the ac- 
tion, these are the well-known Regge equations 

/ 

a discrete version of the continuum Einstein equations in 
vacuum, = 0, namely of the vanishing of the Ricci 
tensor. They give a relation between the deficit angles 
of different faces. The integral (|2"5|) is dominated by its 
stationary 'trajectories', namely by the sets of lengths l s 
which are a solution of the Regge equations (|2"9"|) . However, 



in order to pick up a single classical 'trajectory' we need 
to specify appropriately the boundary conditions and pass 
to the transition amplitudes. The spinfoam boundary for- 
malism for the transition amplitudes is briefly reviewed in 
the next section. 

General boundary amplitudes. Given a 2- 
complex with boundary, its boundary graph T is an ab- 
stract oriented graph made of links I (where the external 
faces end) and nodes n (where the external edges end). 
The boundary graph inherits its labeling from the external 
faces and edges. The set formed by the boundary graph T, 
the spins ji associated to the links, and unit vectors n„i as- 
sociated to the nodes is the boundary data. The spinfoam 
transition amplitude for the 2-complex a with boundary 
T in the Bloch coherent state basis [23] is a functional of 
the boundary datg0 defined as 




where now the action is split into a boundary action plus 
a bulk action. This formula needs further clarifications. 
First, the boundary action Si contains the external face 
amplitudes Pi . For an external face, the formula for Pi is 
the same as ([7]) except that for the edges ending on the 
boundary we have only 'half of ([8]). Second, the sum- 
mation is over the sole internal spins jf £ N, and the 
integrals are over the Spin{A) gauge variables and SU{2) 
variables (the unit vectors) labeling the internal edges. 

The semiclassical analysis of the previous section is 
straightforwardly generalized to a 2-complex with bound- 
ary, which allows to select the continuum/semiclassical 
regime in the following way. First, write the boundary 
functional W(a;, n ne ) in terms of the areas. Then we look 
at the behavior when all the boundary areas a/ are macro- 
scopic, ai 3> lp, in the limit 7 — > with ai fixedO] 

This regime is formally the simultaneous continuum 
7 — s- and semiclassical limit H 0. Now the bound- 
ary amplitude gets its dominant contribution from the so- 
lution to the Regge equations (|2"9"|) compatible with the 
specified boundary data (see [37][33] on the initial value 
problem), namely we have the oscillatory behaviour 

W{a h n ne ) ~ e iS »( ai ' n ™\ 7 -> 0, h^O (31) 

where now Sr is the Hamilton function, that is the ac- 
tion evaluated at the classical trajectory determined by 
the boundary data[f] The set of boundary data can be 

6 The amplitude 1 130 II for a 2-complex with boundary has the nice 
interpretation as transition amplitude associated to a spin-network 
supported on the boundary graph T. 

7 We suppose also that the path integral is dominated by areas 
of the same order of magnitude of the boundary areas: a; ~ af. In 
other words, the macroscopic boundary state must enforce macro- 
scopic areas in the bulk of the triangulation. This is the condition of 
validity of the WKB expansion that has to be checked case by case 
for the specific geometry chosen. 

8 If the boundary data arc not consistent with a boundary trian- 
gulation, the amplitude is exponentially suppressed. 
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equivalently mapped to the set of lengths of the boundary 
triangulation. This set of lengths is a Dirichelet boundary 
condition for the Regge equations of motion to be used to 
determine the classical solution in the interioid and eval- 
uate the Hamilton function in (|31|) . 

The result pip is consistent with what we expected from 
a theory of quantum gravity, and is the concrete realiza- 
tion of the equation ^) in the introduction. 

Conclusions and outlook. — It this letter we have 
discussed a proposal for the semiclassical limit of spin- 
foams truncated to an arbitrary, finite triangulation (2- 
complcx), where most calculations are done. We find 
(equation ([3"TjO that the transition amplitudes are pro- 
portional to the exponential of the Hamilton function of 
Rcgge-Einstcin general relativity, as expected, up to h- 
corrections and 7-corrections in the simultaneous semi- 
classical and continuum limit. The first corrections corre- 
spond to the standard WKB expansion of the path inte- 
gral. The latter are new and are the effect of the discrete- 
ness of geometry, in the sense that the spectra of areas and 
volumes are discrete and the discreteness is controlled by 
the Immirzi parameter 7. The continuum limit we take is 
a looser concept of the 'full' continuum limit defined as the 
phase transition to a smooth spacetime manifold. How- 
ever, working with a fixed, finite triangulation the only 
way of taking a continuum limit is to look at the contin- 
uous spectrum limit of the fundamental geometric opera- 
tors. We have done this by letting 7 run to zero, keeping 
fixed the macroscopic areas a/ (this is still a tentative 
proposal and must be further investigated). Remarkably, 
as explained in the previous section, the result sheds new 
light on the previous calculation of the graviton propa- 
gator with the 'new models' [T5], where the same kind of 
7-corrections to the standard pcrturbative tree-level prop- 
agator has been found (one could speculate on potentially 
observable signatures of those pure LQG corrections). 

It is also very interesting to notice that essentially the 
same continuum limit was considered by Bojowald |36j in 
the context of loop quantum cosmology. Quoting its ab- 
stract: "standard quantum cosmology is shown to be the 
simultaneous limit 7 — > 0, j — > 00 of loop quantum cos- 
mology^ , a strong analogy it is worth studying further. 

We have disregarded other possible contributions to the 
amplitude (symmetry related spacctimcs, vector geome- 
tries, degenerate geometries etc.) that could spoil the cor- 
rect semiclassical behaviour. For example we expect an- 
other term in the leading asymptotics (|31[) of the transition 
amplitudes which corresponds to a spacetime with oppo- 
site orientation that would change the (oscillatory part of 
the) amplitude into the sum of two sign-reversed exponen- 
tials, namely 



as in the asymptotic formula for the amplitude of a sin- 
gle vertex [2"5ll3T)rt4"l"] . In principle this does not pose a 
real problem: for example, see E2] on the fate of the sign- 
reversed exponentials in the bulk of the triangulation, and 
how a coherent boundary state peaked on the appropri- 
ate extrinsic curvature is able to select a single exponen- 
tial (see also E3]). Finally, we have also disregarded the 
potential divergencies associated to bubbles in the foam 
for which a suitable regularization and renormalization 
scheme [44H47] could be required. More details, including 
the Lorentzian signature, can be found in a longer version 
SHI of this letter. 
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